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I . Introduction 

The goal of the first year effort was to calculate scattering 
from cui inhomogeneous layer with irregular boundaries to model 
natural terrains such as a layer of vegetation or sea ice. The 
inhomogeneities were modeled by spherical or disc-shaped discrete 
scatterers which were small compared with the incident wavelength 
and were in the far field of one another. It was found that the 
cross -polarized scattering was dominated by multiple scattering 
effects and was sensitive to the orientations and distributions of 
the scatterers. This model has been applied to interpret measure- 
ments from vegetation, snow and sea ice. 

Tiie goal of the second year was to extend the scattering model 
developed in the first year to handle disc-shaped scatterers which 
are comparable to the incident wavelength and to use the scattering 
model to investigate the relative merits between active versus 
passive sensing of soil moisture over vegetated terrain. Results 
indicate that scattering measurements are more sensitive to soil 
moisture changes than emission measurements. This is because 
wnile both types of measurements lose sensitivity to soil moisture 
because of the vegetation layer, the loss is greater for passive 
than active measurements 


II. Work Accomplished 

The results of the investigation during the first year have 
been documented and published. Details are given in three papers 
appended to this report as Appendices A, B, and C. 

Part of the results of the second year effort has also been 
published. This part appears in Appendix D. Other results of the 
second year study have been documented for publication. It is 
appended as Appendix £ in this report. 


III. Conclusions 

Significant progresses have been made in the previous years. 
We are looking forward to continue our efforts at the University 
of Texas at Arlington to make further progress in the theoritical 
modeling area. 
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Abstract 

A general theory of intensity scattering from small 
particles of arbitrary shape has been developed based 
on the radiative transfer theory. Upon permitting the 
particles to orient in accordance with any prescribed 
distribution, scattering models can be derived, fty 
making an appropriate choice of the particle size, the 
scattering model may be used to estimate scattering from 
media such as snow, vegetation and sea ice. For the 
purpose of illustration only comparisons with measure- 
ments from a vegetated medium are shown. The difference 
In scattering between elliptic- and circular-shaped 
leaves is demonstrated. In the lpw-f requenev limit, the 
major factors on backscatterlng from vegetation are 
found to be the depth of the vegetation layer and the 
orientation distribution of the leaves. The shape of 
the leaf is of secondary importance. 

I. Introduction 

The scattering from randomly or 'anted circular 
discs has been used to model weve scattering from a 
vegetation canopy. 1 * 5 

This study deals with scattering from randomly 
oriented scatterers of arbitrary shape in the low-fre- 
quency limit. First, a general representation of tha 
polarizability tensor in the reference frame is obtained 
in terms of the angle* of orientation and the polariza- 
bility tensor In the local frame (the principal frame 
of the scatterer). Then the polarizabi I i ty tensor is 
used to calculate the scattering amplitude when the 
scatterer Is illuminated by an incident plane wave. 

From the scattering amplitude the backscatterlng cross- 
section for a layer of randomly oriented scatterers may 
be computed based on the first-order solution of the 
equation of transfer. 6 As an i I lustration, the study 
Is specialized to the case of randomly oriented elllp- 
t leal 1 y shaped leaves to model scattering by a vegeta- 
tion canopy. Differences in scattering between circular 
versus elliptic leaves are illustrated. 

2. The Polarizability Tensor In The Reference Frame 

For a wave incident upon an arbitrarily shaped 
scatterer with all dimensions small compared with the 
electromagnetic wavelength, the scattered field due to 
the presence of the scatterer can be represented by that 
radiated from an equivalent dipole with moment F, 7 


where t\ ts the applied field in the absence of the 
scatterer and 5 is the polarizability tensor. 

For a scatterer with its principal axes coinciding 
with the local frame, the polarizabi I i ty tensor can be 
written as, 


, ii x i x i 


where a}] 1 * depend on the electricproperties^and the 
dimensions^ the scatterer,® and xl*s (with xj*x* , 
xj-y* , x>z') are^the local frame axes related to the 
reference frames Xf through angles of rotation (a,6,y) 5 


where Xj m »s are given in Appendix A of Ref. 9. 

To transform the polarizability tensor from the 
local frame to the reference frame we multiply Eq, 2 
from both sides by the unit dyadic 


T £ * " 

' " v. 


yi.ldlng 


■ t i •, 


mn x m x n 


l X lo X !n 


From Eq. 6 it is clear that 5 p is the symmetric tensor* 
The explicit values of ^ mn ’s are given in Appendix A of 
Ref. 9« In the special case of a sphere 7 we have 

a 3(e?-I) 

°i [ ^ ®« ■ — ? v o (7) 

11 4 (eS2) 0 

where e r *e{.+je , f J is the relative dielectric constant of 
the scatterer with respect to that of the host medium; 
and Vq is the volume of the scatterer. (Note: there 

Is a Air difference between values of <*};'* and those 
reported in Refs. 7 and 10. This factor depends on how 
we express the scattered field in terms of ajj's.) 
Substituting Eq. 6 into 5 we obtain 


Vln " V 


where is the Kronecker delta. 

3. The Scattering Amplitudes In The Reference Frame 

For a wave incident upon a randomly orientfd scat- 
terer in i direction with polarization vectors h* and 
v* (Fig. 1), the scattering amplitude matrix for the 
scattered wave in s direction with polarization vectors 
h $ and v s , Is given by 1 




E 5V 

q-v,h 


t. X lm x 


* Z Z fp.tf.oW (9) 

p«v,h q»v,h 

where f p q(s,?) is the scattering amplitude for the 
scattered wave in s direction with polarization p due 
to an Incident wave In ? direction with polarization q 
(p«q«v or h) ; and k is the wave number of the host me- 
dium. The polarization vectors v ! and fi* are related 
to the incident direction I by 

h 1 ■ zXi/|zXi| A £ h|x m - (ycosp.-xsin*.) 


cost), (xcosb I + ys in$j)-zsinb. 





Figure 1. The polarization vectors for tha 
incident and scattered wave. 
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wave 


Similar relations can be written for the scattered 
Substituting Eq. 5 into 9 we can write f p q(i,?) 


I) 
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To show that^Eq. 12 obeys the reciprocity theorem we 
interchange s and ) in Eq. 12, yielding 


V 1 **) 


w 2 3 
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It is worth noting that depends on a, 0, y and c^, n 
but is independent of s and ?. 


k. The Extinction Coefficient For 
Randomly Oriented Scat te re rs 


For a plane wave propagating in i direction defined 
by &((${,$() through a collection of randomly oriented 
scatterers, the total toss per unit length is given by 
the extinction coefficients as, 


s ,! > ■ V ! > <"> 

The subscript a stands for absorption, the subscript s 
for scattering and the subscript q»v or h represents 
the polarization of the incident wave. 

The scattering coefficients can be written as 2 


V'' ' p£h L *• 


(15) 


where dft s is the differential solid angle in the direc~ 
tlon of the scattered wave s, and the ensemble average 
<> is taken over the spatial distribution and orienta- 
tions of the scatterer. The spatial averaging is equi- 
valent to multiplying Eq. 1$ by the number density Ng. 
From Eq. 12 we can write: 




i t £ £ 

ro*J n*l £* I j*l 


1 nw a tJ*M P m P * 

06) 


where * Is the complex conjugate symbol. From Eq. 10 
and Eqs. (A- 1 1)- (A- 17) of Ref. 9 it follows that: 
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Substituting Eqs. 17 and 18 in 16 and then 16 into 
15 we get 

1 2 


k sq (,) 
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cn 
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wh«r« |o^| 2 Is th. value of |i m( .| 2 after Integration with 
reipact to a over (0,2r). After Integrating Eq. 19 with 
raspact to dft s and summing ovar n and scatttrad polari- 
zation we can write the scattering coefficient as: 


N o k / h 
r‘(r 
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Substituting Eqs. (A- 1 2) - (A- 1 7) Into (A- 1 9) of Raf. 9 
and than substituting the result (with l*m and j«n) Into 
Eqs. 20 and 21 we can write the scattering coefficient 
in terms of B, y and a) f as: 

.2.,. , ,2 
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(22) 


k #h (i) » <cos 2 S|a|j i 2 +(sIn 2 B+cos 2 Ycos 2 0)|a^ 2 | 2 

♦(sin 2 0+sin 2 YCos 2 3) |a^j| 2 > (23) 

In Eqs. 22 end 23 Che ensemble average <> is taken over 
the scatterer orientation with respect to the angles 3 
and y. Under the assumption that al.-aU* " nd ym °* 

22 and 23 reduce to Eqs. (68) and (&9) in Ref. 2. 

The absorption coefficient, for a unit jneident 
plana wave in i direction with polarization q*, can be 
written as 11 

*.q(0 • V C r 4 L 'K*‘ 5 J l)dV> 

v o 

- M 0 kc 'r <^q'S ,)> V 0 {zk) 

where is the electric field inside the scatterer, 
which Is related to the induced dipole moment pq through 
the relation 7 


From Eqs. I and 2 we can write Pq is 
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Upon substituting Eq. 26 into 2 5, th« field Inside the 
scatterer Is 
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From Eqs. 24 and 27 tha absorption coafflciant It 
k iq (f) " N 0 k 


In Eqs, 31 end 32 It end -I are In (Oj.nm) and (0p$|) 
directions, rtsptct i ve ly ; d Is the d«pth^of the layer 
containing the scattertrs; r p (f) and r q (i) are the 
Fresnel reflectivities for the Interface separating the 
layer and the lower half-space (p-q«v or h) * From Eqs. 
10, 16, 16 and Eq. (A-19) of Ref. 9 we can write 

<|f vv (-U)i 2 >-N 0 ji <cos\|^| 2 


( 28 ) 
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where 


+sln**0 1 1 a^| 2 > 


k » 


kc*; 


* V 0 |e r -I| z 

Substituting Eq* 10 and Eq. (A-2) of Ref. 9 Into Eq. 28 
and integrating with respect to a over (0,2tt) we obtain, 
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(33) 
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+|a^ 3 | 2 (sln 2 6cos 2 Y^sin 2 Y) > 


(36) 

(29) The explicit forms of Eqs, 33“3$ In terms of 6 and y es 
well as c^ n are given in Appendix C of Ref. 9. 

For an elliptic scatterer with semi-axes a; (1*1 f 
2, 3) the potarizabf I Ity tensor elements a|j arc 7 


(30) where 


For an elliptic scatterer in the low-frequency limit 
Eqs. 29 and 30 are equal to the extinction coefficient 
calculated by using the forward scattering theorem. 0 

5* The Backscatterf nq Cross-Section For A 
Layer 0^ Random! yOrien ted Ellipsoids 


and 
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For a collection of identical sparsely distributed . - 

randomly oriented lossy dielectric ellipsoids embedded ,n Eq * 37 V 0 i# tht vo,umt of tht icatt#rtr 9 |ven b V 

in a layer without upper boundary above a half-space 
Me like- and cross-polarized backscattering cross- 
sections per unit area computed from the first-order 
solution of the radiative transfer equation are 6 
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The integral Aj in Eq. 38 can be found in Ref. 12. In 
case of a thin elliptic disc (ajMjiej) which can be 
used to model a leaf, we can write*2 

A, . liVTT 2 P("/2»«)-E( ft ^.«) 
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wher. .Wl*(«2/*l) 2 • F(*/2,e) and E(*/2,e) ar* complete 
elliptic integrals of the first and second kind, respec- 
tively, 13 given by 
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Figure 2. Backscattering from a 
half-space of randomly orlanted 
elliptic discs [f-l. I GHz, 
c r -30.8+jl.8, ai*2.8 cm, a 2 *0.84 cm, 
•3-0.375 m 9 0*<6<30 # , 0 # <y<90\ 
0*<a<360 # J. 


Figure 3* Sackscattering from a 
half-space of randomly oriented 
el I Iptic discs [f«1 . 1 GHz, 
c r -30.8+J 1 .8, ai-2.8 cm, a 2 *0.84 cm, 
a3-0.375 mm, 60*<e<90 # , 0 # <y<90’, 
0 j <a<360«]. 


Figure 4, Backscatter ing from a 
half-space of randomly oriented 
el I i; tic discs [f-l .1 GHz, 

C--30.8+J I ,8, a»*2.8 cm, a 2 *0,84 cm, 
•3*0,375 mm, 0 # <B<90\ 0*<y<30 # , 
0*<ot<360 # J. 



Figure 5. Backscattering from a 
ha If- space of randomly oriented 
elliptic discs [f-l. I GHz, 
cr" > 30.8+jt*8, at-2.8 cm, 02*0.84 cm, 
•3-0.375 mm, 0 # <8<90*, 60*<y<90\ 
0Va<360*J. 


and 


fit/ 2 


V^7 


In 2 S 


dS 


Figure 6. Backscattering from a 
half-space of randomly oriented 
elliptic discs [f-t.l GHz, 
c r -30.8+jl.8, aj-2.8 cm, a 2 *0.84 cm 
•3*0.375 mm, 0*<8<90* , 0 # <y<9C # , 
0^<a<360*J. 


Figure 7* Cross-polarized pattern 
at normal incidence along circular 
disc axis [f-l. I GHz, c r *30.8+j 1 .8, 
•j-a 2 *l cm, 03-0.872 mm, 0*<a<j60*J. 


E (ir/2, e) - / 

'0 

When 0|*a 2 (oblate spheroid), we can write 22 
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where ro-ai/aj. 
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aj, *2* * nd 03 are taken to be parol - 
axes, respectively. This means that 
the orientations of the semi-axes •], a 2 » 03 are de- 
scribed by (a,6) , (a,8,y) and (o,8,y). 9 The tilting of 
•I and a 2 with respect to the local axes are of course 
controlled by 8 and y, respectively. 


6. Comparisons Between Circular And Elliptic Discs 

In Figs. 2-6 backscattering characteristics from 
elliptic discs, with dielectric constant based on for- 
mula given in Kef. 14, are shown for five different com- 
binations of orientations. In ail cases the orientation 
distributions of the axes not specified below are as- 
sumed uniformly distributed over (0*,90*)« 

(1) The major axis Is restricted to small 
angles, i.e., 0*<8<30 9 . 

is restricted to large 


is restricted to small 
is restricted to large 


tilt 

tilt 

tilt 

tilt 


(2) The major axis 
angles, i.e. r 60*<8<90*. 

(3) The minor axis 
angles, i.e., 0 # <y<30*. 

(4) The minor axis 
angles, i.e., 60*<y<90*. 

(5) Both axes are uniformly distributed over 
(0\90 # ). 

Before discussing the characteristics of the cases 
shown, note that polarized scattering is usually in 
agreement with intuition. To understand cross-polarized 
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Figure 8«i Cross-polarized pattern at normal Incidence 
along major axis of an elliptic disc [f»l.l GHz, 
c r »30.8+jt .8, at«2.8 cm, a2»0.84 cm, a^O.375 mm, 
0 V <ai360 # ] . 

scattering It should be remembered that the cross-polar- 
ized reradiation pattern for small circular discs In- 
creases with the increase In the local incidence angle 
(Fig. 7); while for an elliptic disc of the same volume, 
this reradlation pattern could have a dip with the in- 
crease In the local Incidence angle along the plane 
parallel to the major axis (Fig. 8a). The reradlation 
pattern parallel to the minor axis shows a minimum at 
near normal incidence (Fig. 8b). 

For case (I) (Fig. 2) the horizontally polarized 
component at large incidence angles is higher than the 
vertically polarized component due to the more nearly 
horizontal disc orientations. The reverse is true in 
case (2) (Fig. 3) as expected. For the cross-polarized 
component case (1) shows a slower angular drop off than 



Figure 8b. Cross-polar! zed pattern at normal incidence 
along minor axis of an elliptic disc [f»!.l GHz, 
c r "30.8+jl .8, aj-2,8 cm, a2*0»84 cm, ai*0.375 mm, 

0*<a <360*1 . 

case (2), since the cross pattern along the major axis 
or minor axes is close to its maximum at near grazing 
Incidence and assumes smaller values at near normal 
Incidence. In cases (3) and (4) (Figs, 4 and 5) the 
change In distribution is associated with the minor 
axis. As expected, the polarized reradlation pattern 
is more Isotropic over the minor axis than the major 
axis. Hence, this change in distribution produces only 
small changes in the polarized scattering components. 
For the cross-polarized component, the smaller tilt 
angles [case (3)1 lead to smaller returns near normal 
incidence and a slower angular drop off than case (4). 
This Is expected since the cross pattern is at a mini- 
mum near the local normal of the disc parallel to the 
minor axis (Fig. 8b). These are only small differences 
between cases (4) and (5) for the polarized component 
for the same reason given for cases (3) and (4). The 
cross-polarized return of case (5) lias between those 
of cases (3) and (4) as expected. 

In Figs. 9-11, the backscattering curves for cir- 
cular discs of the same volume as the elliptic discs 
with major axis equal to the minor axis, the five cases 
discussed in the previous paragraph reduce to three 



Figure 9« Backscattering from a 
half-space of randomly oriented 
circular discs [f«t.l GHz, 
e r »30.8+j I .8, ai*a 2 "l cm. 33»0.872 
nw, 0*<B<30*, 0*<y<90 # , 0*<a<360*] . 



Figure 10. Backscattering from a 
half-space of randomly oriented 
circular discs [f-l.l GHz, 
e r »30.8+jl.8, a | I cm, a*"0.872 
mm, 60*<S<90\ 0*<y<S0\ 0*<a<360*]. 
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Figure II. Backscatter ing from a 
half-space of randomly oriented 
circular discs [f»l.t GHz, 
e r »30.8+jl . 6 , ai»a 2 *l cm, aj»0.872 
mm, 0*<B<90*, 0*<v<90*, 0*<o<360*]. 
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Figure 12. Comparison between 
beckscattering *rom a half-space, 
a layar with plana and rough 
bottom boundarlas [f*I.I GHz, 
c r *30.8+jl.8, c fl *IO, ka» I • I » 
kt*8 , VaNo"0.002, 60 # <8<90\ 
0 # <y< 90*, 0*<a<360*J "* 



Flgurt 13* Comparison with measure- 
mants of whaat (f*l.5 GHz, 
c r *24.5+JI»7, c-*2$, o«l.l cm, 1*25 
cm, a|*a2"l»3 cm, 03*0.02 cm, 
0*<8<90\ 0*^y< 90 # , 0 # <a<j60*]. 



Figure Ik. Comparison with measure- 
ments of whaat [fnk.2$ GHz, 
e r *23.5+j4.7, c-*25, o*l,l cm, t*IJ 
cm, a|"a2*l.3 cm, 03*0.02 cm, 
0*£8<90 # , 0*<y£90\ 0 c <c£360 # ], 


casas. Figs. 9, 10 and II ara comparabla to F tgs. 2 or 
4, 3 or 5 and 6. Tha polarlzad and cross-pot (zad 
scattaring also follow tha ganaral pattarns In tha cor- 
responding figures for reasons similar to those given 
in tha previous paragraph. Of course, small differences 
exist in level and trend between the corresponding 
figures especially for the cross-polarized return (see 
for example, Figs. 2, 4 and 9)* In view of tha dif- 
ferences between Figs. 7 and 8a, one would expect much 
greater difference between elliptic and circular discs. 
However, due to the assumed random orientation, the 
difference due to disc shape in this low-frequency case 
Is not as pronounced as one might expect. Instead, 
significant differences In scattering due to different 
distributions are very obvious in the cross-polarized 
return. This Implies that the cross-polarized return 
may serve as one contributing factor for crop identifi- 
cation when significant differences exist between the 
orientation of leaves and when the leaf distributions 
In orientation are known for vegetation types to be 
different iated. 


fairly flat for the cross-polarized return and peaked 
near nadir for the polarized return indicating the 
strong influence of the rough ground Interface, Al- 
though the leaf of wheat Is elliptic, use of circular- 
shaped leaves with random orientations appears to pro- 
duce good agreements with data. This means that the 
effect of leaf shape is Indeed not an important factor 
In low-frequency scattering. 

8. Conclusions 

In vegetation modeling at low frequencies, the im- 
portant factors governing the backscattering curves are 

(1) the roughness of the vegetal ion-ground interface, 

(2) the leaf orientation distribution, and (3) the depth 
of the vegetation layer. It is found that the effect of 
the leaf shape Is not as Important a factor although It 
does have some effect on the level and trend of the 
angular backscattering curve. 
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7. layer Effect And Comparison With Measurements 

In Fig. 12 comparisons are shown between backscat- 
*.erfng curves from a vegetated half-space, a vegetated 
layer without boundaries, a vegetated layer with plane 
bottom boundary and a vegetated layer with rough bottom 
boundary. It Is seen that scattering from a half-space 
has a significantly slower angular trend than that of a 
half-space. For horizontal polarization, the difference 
between a plane versus a rough layer boundary Is small 
at large angles of Incidence for polarized scattering 
and at smalt angles of incidence for cross-polarized 
scattering. 

Since the effects dominating the backscattering 
curves at low frequencies (scatterer small compared with 
the electromagnetic wavelength) are (I) leaf distribu- 
tion, (2) layer versus half-space, and (3) surface 
roughness, In Figs. 13 and 14 we show comparisons with 
backscatter measurements from wheat 15 using a vegetated 
layer of circular leaves with a rough bottom boundary. 

The rough surface scatter model used is the Kirchhoff 
model under scalar approximat ion 26 characterized by a 
standard deviation of surface height o and surface cor- 
relation length i. In both figures, the data trends are. 
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A vegetation layer is modeled by a collection of randomly ori- 
ented circular discs over a half space. The backscattering coef- 
ficient from such a half space is computed using the radiative 
transfer theory. It Is shown that significantly different results 
are obtained from this theory as compared with some earlier inves- 
tigations using the same modeling approach but with restricted disc 
orientations. In particular, the backscattered cross-polarized 
returns cannot have a fast increasing angular trend which is incon- 
sistent with measurements. By setting the appropriate angle of 
orientation to zero the theory reduces to previously published 
results. Comparisons are shown with measurements taken from mllo, 
corn and wheat and good agreements are obtained for both polarized 
and cross-polarized returns. 


1. INTRODUCTION 

The problem of scattering from randomly 
oriented scatterers with application to 
vegetation has received mush attention in 
recent years. The radiative transfer 
theory has been used to calculate the back- 
scattering from a layer of randomly ori- 
ented ellipsoids with application to vege- 
tation [Tsang et al. , 1981]. Foldy's 
method and the distorted Born approximation 
have been used to calculate the backscat- 
tering cross section for randomly oriented 
circular discs in a half space [Lang, 1981]. 
The major advantage of this approach is 
that most of the model parameters are di- 
rectly related to ground truth measurements. 

In this study we use the scattering am- 
plitude formula by Karam and Fung [1982a] 
for randomly oriented scatterers of arbi- 
trary shape to derive the scattering am- 
plitudes for randomly oriented circular 
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discs. Then the phase matrix elements and 
the extinction coefficients are derived. 
These quantities are needed in the first- 
order radiative transfer theory for com- 
puting scattering from a half-space medium 
containing randomly oriented sparsely dis- 
tributed lossy dielectric discs used to 
model a vegetated medium. In this model 
the effect of air-vegetation boundary is 
neglected due to the sparse distribution 
of the leaves. The decay of the incident 
wave in the vegetated volume is accounted 
for by the extinction coefficient which is 
related to the scattering amplitudes in 
the forward direction [ishimaru anu Cheung, 
1980; Karam and Fung, 1982b]. 

The main difference between this study 
and those of Tsang et al. [1981] and Lang 
[1981] is that our formulation permits the 
circular discs to be oriented in any direc- 
tion and accounts for changes in polariza- 
tion of the scattered field due to changes 
in orientations of the circular discs. 

This difference is illustrated numerically 
in section 5 where comparisons with mea- 
sured data on wheat, corn and milo [Dobson 
et al. , 1977] are also shown. 
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2. THE SCATTERING AMPLITUDE MATRIX AS * - cos6 k - sin9 v (5c) 

A SPECIAL CASE OF THAT OF ELLIPSOIDS X 1 


Consider « randomly oriented ellipsoid 
with its semlaxc^ (a,b,c) aligned along the 
coordinates (x'.y’.z') of a local frame 
related to the principal frame (x,y,z) 
through the angles of rotations^ (a,8,y) 
with respect to the z', y ’ and x' axes, 
respectively (called the xyz convention by 
Goldstein (1981]) [Karam and Fung, 1982a], 
Thus, 

x' • cosScosax + cosBsinay - sinez (?) 

y* ■ (cosasinBsiny - sinacosy)x 
+ (cosacosy + sinasin3siny)y 
+ cosBsinyz (2) 

■ (sinasiny + cosasinScosy)x 
- (cosasiny - sinasin6cosy)y 
+ cosScosyz (3) 


Then substituting (5a)-(5c) into (3) we get 

*’ ■ 'kA + ‘vi’l + 'hi 6 ! <6) 

where 

c ki " (sin(o-4» i )siny + coste-^sinBcosy] 

* sinQ^ + cosBcosycosQ ^ (7a) 

t yi ■ (sin(a-^ i )siny + cosfa-^sinBcosy] 

• <5036^ - cosBcosysine^ (7b) 


t hi * sin(a-4> i )sin0cosy - cosU-^siny 


Now substituting (6) into (4b) we get 
* , 1 


1 n — 

VSi + 


2 

hi 


! Vi - c hi h i ] 


(7c) 


(8a) 


For an incident plane wave in the k^ di- 
rection the polarization vectors (v^ and 
h^) in the principal frame, and the cor- 
responding polarization vectors in the 
local frame (v^' and h^’) can be expressed 
as 

h A - z X itj/lz / k ± | v A ■ h 4 X k t (4a) 

h i ' - i* x k 4 /|i» x kj v t » - h i * x k ± 

(4b) 

where 

k^ ■ sinSjCosSjX + sinfl^sin^^y + cosB^z 

(4c) 

To relate h^', v^' to hj, v^^we express^the 
principal frame in terms of k^, v 4 and 
as 

A A A A 

x ■ sine^cos*}!^ + cosBjCosiji.jV^ - sin^h^ 

(5a) 

y ■ slnO^sln^k^ + cosQjSin^v^ + cos^h^ 

(5b) 


V • 




2 2 
t + t 

C vi * c hi 


>Vi 


- V h tl 


(8b) 


Relations similar to (4)-(8) can be written 
for the scattered wave by replacing the 
subscript 1 with s. (Note that the ori- 
entation of the scatterer and its effect 
on the polarization of the scattered field 
are not represented in the most general 
form when the local frame is described 
by the Eulerlan angles of rotation (appen- 
dix) . ) 

Under the low frequency approximation 
[Stratton, 1941; Stevenson, 1953; Ruck et 
al., 1970; Ishimaru, 1978; Tsang et al., 
1981] the scattering amplitudes for the 
ellipsoid in the local frame can be writ- 
ten as 


f' 

w 


(s,i) 


4 H Vv e > 



1 •x’) (v^' 'X*) 

vA + 1 
a 



>(v '•y*) (v 


v s ’»z')(v i ‘.z')* 

vA„ + i 

C M 


(9) 
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hv 


<••*> * TnfVV C >| 

. v *. * 1 

(v 1 , *y , )(h 1 '*y , > 


v S + 1 + 

vA_ + 1 
c 

- 4Vv c > 

’(h # , *x , )(v 1 '*x') 

vA g + 1 

Ch — * • y* > (v t * ’y') 

(h # ' *z') (v^' *z' ) 

vAb + 1 + 

vA„ + 1 
c 

t (s,i) - %rV 0 (c 8 -e) 


L VA . + 1 

(h ••y , )(h/-y') 

(h 8 , -i , )(fi i , ** , > 

vA. + 1 
D 

vA + 1 
c 


( 12 ) 

where the third terrain (10)-(12) is zero 
since (h^* • z' )"(h s ' ‘z'XO; tu is the angular 
frequency of the incident wave; e 8 is the 
permittivity of the scatterer; c is the 
permittivity of the surrounding medium; and 
Vg"(4Trabc/3) is the volume of the scatter- 
er. In (9) -(12) we also have 

R(u) ■ £(u+a 2 ) (u+b 2 ) (u+c 2 )J 

-■fp 


+u)R(u) 


i a ' 

f 

Jc\ 


du 

+u)R(u) 


a and thickness 2c. In this limit [Tseng 
et el., 1981] we have 


A - A, - A - 

* ^ 2a 3 c a 2 c 


Hence, 

vA . ■ VA b ■ (t - 'll! vA c • (-f - l ) 

Let us define 

(e / e-1) (e / c-1) 2 

« 2 - 2 c - 

*T vA +1 *N vA +1 0 4 71 v 0 

a c 

We can now write the element of the scat- 
tering amplitude matrix for a circular disc 
in the local frame as a special case of 
(9)-(12). 

- C 0 l( V*V )a T 

+ (v s , *z')(v i '"z')(a N -a T )] (13) 

f vh (a * i) " C 0 [( V*V )a T 

+ (v s '*z')(h 1 '*z , )(a N -a T )] (14) 

*!*<■»*> • C 0 [(i ‘s , *V )a T 

+ (h g '*z , )(v i '*z , )(a N -a T )] (15) 


£ hh < * >1) 


c 0 Kh s ' 


‘i'K 


+ (h g , *z')(h i ’*z , )(a N -a T )] 


(16) 


du 

(c 2 +u)R(u) 


When a*b and c<<a the ellipsoid can be 
approximated as a circular disc with radius 


In view of (8) an incident field ampli- 
tude in the principal frame can be expres- 
sed in the local frame where scattering 
takes place in accordance with (13)— (16) . 
Then a similar relation to (8) may be 
applied to convert the scattering field 
back to the principal frame [Karam and 
Fung, 1982a]. Thus, 


f w (8 » i) " Ks f w <s * i) S,i- t vs f vh <s ' 1)t 


hi 


“ t hs f hv ( ®’ l) t vi +t vs f hh (a ‘ C hi ] /D( ® ’ l) 


( 17 ) 


<1 

4 


r;.«. A- 

OF POOR QUALITY, 
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V1*W V * ' "hi T "vs‘vh' "vi 


" t h« f hv (,,£) t hi“ t h» f hh ( * ’ ^ f v J /D(s » ^ 


f hv (# » 1) " It h. f w (l ' i)t vi +t vs f hv (8 ' 1)t vl 
’ c h» £ vh ( ^ c hl“ C vs f hh ^ ^ e hi ^ /D ^ 


£ h h (t,i) m It h- f i*( 8 ' i ) t h< +t v. f h«<»» 1 ) t 


hs" w va ’ "hi r "vs‘hv'“ ’ "hi 


+t hs f vh (s *^ t vi +C vs f hh (s,£)t vi^ /D(s,£) 


where 

D(s,i) 


[mtti+t.)] 


EquaCions (17)-(21) will be used in Che 
following sections Co calculate the first- 
order scattering cross section for a half 
space of randomly oriented circular discs. 

3. FIRST-ORDER SCATTERING CROSS SECTION 
FOR A HALF SPACE OF RANDOMLY ORIENTED 
CIRCULAR DISCS 


For a collection of identical randomly 
oriented lossy dielectric circular discs 
embedded in the lower half space, the 
backscattering cross section [Tseng and 
Kong, 1978; Tsang et al. , 1981] based on 
the first-order solution of the radiative 
transfer equation [Karam and Fung, 1982c] 
can be written as 


a a 9 tr 

’T 


cosd 


A A 


w' 


w' 




cose 


Im<f hh (i,i)>+Im<f hh (-i,-i)> 


(18) 




(24) 


(19) 


(20) 


( 21 ) 


where X is the wavelength of the Incident 
wave, Im( ) is the imaginary part operator 
and the angle brackets are the ensemble 
average symbol over the orientation and 
distribution of the scatterer. 

4. THE ENSEMBLE AVERAGE OVER THE ANGLES 
OF ORIENTATION FOR FORWARD 
AND BACKSCATTERING 

A In the forward direction* k s «k^; hence, 
h s "hi» v s" v i* v s '-v^' and h s '*hi'. Similar 
statements can be made about t V i, t^i and 

c vs * c hs * 

In the backscattering direction similar 
results are obtained as in the forward 
direction except that h s ' and tj, s will dif- 
fer in sign from those in the forward di- 
rection. 

From the above discussion the forward and 
backscattering scattering amplitudes in the 
local frame (13)- (16) reduce to 

■ C 0 [a T +(C 1 ’-;') 2 ( V a I )] (25) 

(26) 

• lc o[v< 6 i'- i ') J, V*i>] < 27) 


f ;„(=!.!) ■ - 0 


Also, the scattering amplitudes in the 
principal frame described by (17)-(21) 
reduce to 


< |f vv (*i.i)| > 


, j 2 tt 

°vh < - 1 * 1) "T 


cos6 . 


Im<f vv (i,i)>+Im<f hh (-i,-i)> 


f w (±i f) - C 0 [v4<V a T>] 

(28) 

i)> f vh (± * ,£) " ^fchiSri^V] 

(29) 

<22) f hv (± *’ £) " ±C o[ t hi t vi (a N" a T > ] 

(30) 

f hh (±£,£) " C o[ a T ft hi (a N” a r ) ] 

(31) 


< l f vh ( " 1>i) l ' 


(23) 


To calculate the backscattering cross 
section in (22)-(24) we first substitute 
(7) in (28)- (31) and then we calculate 
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the ensemble average of Che quantities 
<|f hh <-l.i)l 2> . *n<i <|f vh <-l,l)l 2 >. 

Due to Che symmetry of Che dlec we will 
assume an equally likely diacribuclon func- 
Clon wich reapecc Co a. The ensemble aver- 
age over a will remove che $£ dependence, 
yielding 

,w/2 

P < 8) dS 


<f vv (ii,±i)> - C Q n 0 


H 


sin^cos^j+cos^sin 2 © 


P(Y) *{($ 

■ cos 2 Y+|cos 2 e i sin 2 Y|(a N -a, r ) 


i) 


(32) 


< |£ hh <-i i i)| 2 >.C 0 2 


r r n 
n o[l* T l 2+ J i 


p(B) dS 


2 2 2 
ain Scot Y+sin y 


) 


Ml* / * 

•I p(y) dyj^i 

Jo 

R« ^a T *(a N -a T ))+| ^»in 4 Scos 4 Y 
+ ■|ein 2 Ssin 2 2Y+«in 4 Y) | a^l 2 j] 

: o 2n off /2 

~Un 1 


< w-i.i>r> 


p(S) dB 


(35) 


hh' 


r 

/•ff/2 

/•n/2 

9 t / ^ A _ k 

’ C 0 n oK 

.+ I p(6) d$ 

• *1 p(y) dY |(cos 2 6 sin B+sin 2 20sin 2 0 . ) 

L 

Jo 

•'O 


( n ft 

p(Y) dYj-|/sin 2 0cos 2 Y+sin 2 Yj| 

° I 

(a N -a T*J 

r r ff/2 

n o[|a T l 2+ J * 


(33) 


<|f w (-i,i)| 2 > - c Q 2 


p(B) dB 


1 n/2 

' 


P(Y) dY 


in 2 6co* 2 © jj+cos 2 Bain 2 © ^ 

Y+jcos 2 © ^a in 2 y j Re ^a T *(a N -a T )j 
+ |'|co3 4 0 i ain 4 Y+|*ln 4 0 1 coa 4 8 
+ ■|co8 4 6^ain 4 8+^8in 2 20^sin 2 2Bj cos 4 y 
+ ~(co8 4 0 1 8ln 2 B+8in 2 20 1 co8 2 s)8in 2 2Yj 

• Iv>tI 2 |] 


cos 


(34) 


4 2 4 /1 2 2 

• cos y+cos 0^sin Y+I^sin Bcos 0^ 

+ cos 2 Ssin 2 6 1 |sin 2 2Yj|a N -a T | 2 l (36) 


where ng is che number of che scaCCerers 
per unic volume; p(3) and p(y) are Che 
orienCaClon probabilities wich respccC Co 
B and y, respectively: Re( ) is Che real 
pare operacor and che ascerlsk is Che com- 
plex conjugace symbol. The lncegratlon 
wich respecc co B and Y is caken over ir/2 
due co Che symmeery of che disc. From 
(32)- (36) we noce that ng appears in both 
the numerator and denominator of (22)- (24) 
so it will be canceled out. 


5. NUMERICAL RESULTS AND DISCUSSION 


To compare our formulation wich Chat in 
the literature [Lang, 1981; Tsang et al. , 
1981] and to show che agreement between 
this Cheory and some measured data [Dobson 
et al. , 1977] we will assume an orienta- 
tion distribution funedon wich respecc co 
B and y as 
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Fig, 1. A comparison with Tsang at al. 
[1981] (£-1.1 GHz; a-10 mm; c-0.375 mm; 
e r -30.8+jl.8). 


P(G) 


n 2" n l 


n i <JJ<fi 2 


otherwise 


(37) 


with 0-8 and y. Also, we assume that the 
discs have a dielectric constant based on 
the formula given by Fung and Ulaby [1978] . 

In our formulation we have the freedom 
to select different distribution functions 
with respect to 8 and y and we note the 
following: 

1. When we set y equal to zero our for- 
mulation reduces to those derived by Tsang 
et al. [1981] and Lang [1981 ] . In Figure 1 
a case shown in the work by Tsang et al. 
[1981] is reproduced. It is seen that when 
■y is not regtricted^to zero, substantially 
different oyy and o^y are obtained. Since 
naturally occurring leaves are, in general, 
not restricted in their orientation, we 


believe that both Y and 8 should be per- 
mitted to vary. In Figure 2 another case 
is illustrated wglch shows that the rela- 
tive Ravels of oyy and 0 ^ ure interchanged 
and o^h has a higher level and a different 
angular behavior as compared to those re- 
ported by Lang [1981] , when both y and 8 
are assumed equally likely distributed. 

2. The roles of y and 8 may be inter- 
changed without affecting the o° computa- 
tion, if they are assumed to have the same 
distribution (sec Figures 1-5). This is 
expected since 8 measures tilting in the 
x' direction, while y measures tilting in 
the y' direction, respectively. 

3. When the discs are oriented more 



Fig. 2. A comparison with Lang [1981] 
(f-1.1 GHz; a-10 mm; c-0.25 mm; 
e r -31+11.8). 
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nearly horizontal, we get 0 hh > o VV (Figure 
1). When the discs arg oriented aore 
nearly vertical, 0yv >o HH (Figure? 3-5). 

In Figures 3-5, comparisons with measure- 
ments at large angles of Incidence are 
shown (30*-80*). This angular range Is 
selected to avoid the ground effect in 
backseat tering. In the case of corn and 
milo, it is known that these crops have 
large leaves which do bend and twist in 
their natural state. As a result, one 
leaf may have several scattering centers 
and, hence, is modeled by several discs. 

The leaf thickness, however, should cor- 
respond to disc thickness in modeling. Of 
course, the dielectric property of the leaf 
must be directly applicable in modeling 
also. Figures 3-5 show that the leaves of 
these crops are more nearly vertical than 
horizontal and that among the crops con- 
sidered, the leaf volume is the largest for 
corn and the smallest for wheat. Reason- 



Fig. 3. Comparison with measurements of 
corn (f*l.l GHz; a-15 mm; c«0.167 ran; 
e r »30. 8+J 1. 8) . 



Fig. 4. Comparison with measurements of 
milo (f"1.5 GHz; a»12.3 mm; c-0.2 mm; 
e r -28.52+j2.13>. 


able level and trend agreements are ob- 
tained in Figures 3-5 in both like and 
cross polarizations. 

6. CONCLUSIONS 

By permitting arbitrary orientation in 
modeling, reasonable agreements are ob- 
tained between theory and low frequency 
data (1.1-1. 5 GHz) collected from corn, 
milo and wheat. It is shown that when the 
orientation of the leaves is restricted 
[Lang, 1981; Tsang et al., 1981], the 
levels of 0$y and ogj| are interchanged and, 
hence, cannot produce agreement with these 
crop data., 

APPENDIX: THE EULERIAN ANGLES OF 

ROTATION AS A LOCAL FRAME 

The coordinates described by the Eulerian 
angles of rotation are obtained by three 
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Fig. 5. Comparison with measurements of 
wheat (f«1.5 GHz; a*10 mm; c«0.087 mm; 
e r -28.52+J2.13). 


■ ‘ki^i ♦ ‘vl°l + ‘hi 6 ! 


where 

t ki ■ sin6sin6 i cos(a-4i i ) + cosScosQj 
t vi ■ sinScose^costa-^j) - cosBsinB^ 
t hi ■ sinSsinfa-^) 


(A4) 

(A5) 

(A6) 

(A7) 


We see from (A6) and (A7) that they are 
independent of angle y. Hence, when we 
use (8) to represent the polarization vec- 
tors in the local frame in terms of the 
polarization vectors in the principal 
frame, the system depends only on a and 6. 
In fact (A6) and (A7) are a special case 
of (7b) and (7c) when y m 0 and only two of 
the three angles a, 8, and y are indepen- 
dent in the Eulerian description. Using 
the spherical coordinates to describe the 
local frame of the scatterer [Lang, 1981, 
Figure 3] will lead to the same result as 
the Eulerian angles of rotation with r, <P, 
and 6 corresponding to z' , a, and 8, re- 
spectively. 

Another simple way to see the difference 
between the xyz and y conventions is to 
consider the special case 8*0. Then the 
xyz convention in (1)— (3) gives 


successive rotations. around the z axis, the x’ ■ cosax + sinay (A8) 

y axis and then the z axis again with an- 
gles a, 8, and y, respectively (called the y* ■ -sinacosyx + cosacosyy + sinyz (A9) 

y convention by Goldstein [1981]). The 

resulting coordinates can be written as z' ■ sinasinyx - cosaslnyy + cosyz (A10) 

[Arfken, 1970] 

. . and the y convention in (A1)-(A3) gives 

x' ■ (cosycosScosa - slnysino)x 

+ (cosycosesina + sinycosa)? - cosysine2 *' " C ° 8 < a+Y) * + ain(a+y)y (All) 

(Al) y' - -sin(a+y)x + cos(o+y)y (A12) 


y* • -(sinycosScosa + cosysina)x 
+ (-sinycosSsina + cosycosa)y + alnysinSz 

(A2) 

z' ■ slnBcoaax + slnSslnay + cosBz (A3) 


A . A 

z' - z (A13) 

Thus the y convention depends only on one 
angle, a+y, while the xyz convention still 
depends on two angles a and y. 


If we use the coordinates described by 
(Al), (A2) and (A3) as a local frame ’’o 
describe the orientation of the scatterers, 
we can substitute (5) into (A3). This 
gives an expression similar to (6), i.e., 
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ABSTRACT 

A comparison between active and passive sensing of soil moisture over 
vegetated areas Is studied via scattering models. In active sensing three 
contributing terms to radar back scattering can be identified: (1) the ground 

surface scatter term; (2) the volume scatter term representing scattering from 
the vegetation layer; and (3) the surface-volume scatter term accounting for 
scattering from both surface and volume. In emission three sources of con- 
tribution can also be Identified: (l) surface emission, (2) upward volume 

emission from the veget K u,on layer, and (3) downward volume emission scattered 
upward by the ground surface. As ground moisture Increases, terms (1) and (3) 
Increase due to Increase In permittivity In the active case. However, In pas- 
sive sensing, term (1) decreases but term (3) Increases for the same reason. 
This self-compensating effect produces a lo»s sensitivity to change In 
ground moisture. Furthermore, emission from vegetation may be larger than 
that from the ground. Hence, the presence of vegetation layer causes a much 
greater loss of sensitivity to passive than active sensing of soil moisture. 


1.0 INTRODUCTION 


Theories on Intensity scattering [1,2] and emission [3,4] have been 
developed for an Inhomogeneous layer with irregular boundaries. To model a 
leafy vegetation layer above an Irregular ground surface the phase function 
for a dielectric disk Is needed (Appendix A) and the top layer boundary may be 
removed [5]. Such a model Is expected to be valid for vegetated medium where 
scattering Is dominated by leaves. Note that a real leaf Is, In general, not 
a flat disk but Is curved and may twist, especially If It Is long. Hence, a 
dielectric disk will model a scattering center on a leaf and a long leaf may 
have several scattering centers. 

For soil moisture sensing only polarized scattering has been used In 
practice [6]. It Is also known from experimental studies that the albedo of a 
vegetated medium Is usually around 0.3 or less [7,8]. This means that a 
first-order solution of the radiative transfer equation obtained by assuming a 
weak scattering medium can provide useful estimates. The computational proce- 
dure to obtain the first-order solutions for both the active and passive 
problems are outlined In the next two sections. The relevant characteristics 
of the scattering model Is given In .. ctlon 4 and those of the emission model 
In Section 5. Comparisons with measurments are shown In Section 6. 

« 

2.0 THE FIRST-ORDER SOLUTION FOR THE ACTIVE PROBLEM 

Let us assume that the radiative transfer equations for the upward and 
downward Intensities, I + and I", are applicable within the vegetation layer 
(Fig. 1) 
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where K s , K e are the volume scattering and extinction coefficient matrices, 

m 

respectively,; u s * cos0 s , u * cose; and P( ) is the phase matrix. We assume 
that the layer has no upper boundary and the boundary condition at the lower 
boundary is 
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where G is the surface scattering phase matrix given in [2]. To account for 
polarization effects the intensity column matrices contain the four Stokes 
parameters as their elements [9]. Once the total scattered intensity for an a 
polarized component I a s of the intensity matrix at z*0 is found, the 
scattering coefficient for this component is defined relative to the incident 
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The procedure to determine up to the first-order solution of (1) and (2) 
subject, to (3) is as follows: 

(a) Expand P, I into Fourier series with respect to the azimuthal angle 
and consider one Fourier component at a time. 

(b) Convert a Fourier component representation of (1) and (2) into inte- 
gral equations. 

(c) Use an N-point quadrature integration to rewrite the integrals with 
respect to u in (1), (2), and (3) as a matrix product. 

(d) Solve the resulting integral equations iteratively by assuming that 
K $ is small. 

The m** 1 Fourier component of the scattering intensity including only the 
zeroth-order and the significant terms in the first-order is 


I m « [f x + N + + (f x N" + M fj)] I im , 


(5) 


where is the Fourier coefficient of the incident Intensity. In (5), 
fj represents surface scattering attenuated by the layer. 
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where f m *l/2, m*0, and f m *l/4, m>0, and e p 
as 


is the diagonal matrix defined 
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K el * K ev* extinction coefficient for vertical polarization, 

K e2 * K eh* ext ^ n ct1on coefficient for horizontal polarization, 

K e3 * K e4 * < K ei + K e2)/ 2 * 

and p, 1 are the Indices for the Stokes parameters and the quadrature points, 
respectively, and G m is the m th Fourier component of the surface scattering 
phase matrix [2]. ?i + in (5) represents the volume scattering operator which 
scatters a downward propagating intensity upward. Its elements are defined in 
terms of the elements of the phase matrix for the pq polarization component as 


* K spi f ™ t 1 - e' (Ke > 1+K *«i )d ] 


ij 


‘WW 1 • 
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where K spi » K sp /jj if K $p is the volume scattering coefficient for p polari- 
zation, and vij ■ cos0j. The Index j is also for denoting quadrature points 
similar to 1 except^lt Is used for the Incoming polar angle direction. The 
last term In (5), (Mf^fjN'), represents combined surface- volume scattering. 
The matrix elements of M for the pq polarization element is 
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M may be Interpreted as the volume scattering operator which scatters an 
upward propagating Intensity upward. Similarly, N" Is the volume scattering 
operator which scatters a downward propagating Intensity downward. Its 
elements are 
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Since we assume a weak scattering medium, the dominant term In (5) Is 
(6), especially at small angles of Incidence. On the other hand, the surface 
scattering matrix In (6) can be a fast decreasing function of the incidence 
angle e. If observations were made at e>25° It is possible for the other 
terms in (5) to dominate. For soil moisture sensing only small e has been 
used. Hence, the effect of the vegetation layer is mainly attenuation. As 
soil moisture Increases, both the surface and the surface-volume scattering 
terms In (5) increase. Hence, the decrease in the surface scattering term due 
to layer attenuation Is compensated partly by the surface-volume scattering 
terms. 


3.0 THE FIRST-ORDER SOLUTION FOR THE PASSIVE PROBLEM 
* * 

For the emission problem, the radiative transfer equations can be written 

In terms of the upward and downward temperatures T + , T“ and the temperature of 

the layer, T A . To account for polarization, T* is taken as a column matrix 
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I fart 


containing two elements, T v and T^, corresponding to vertical and horizontal 
polarizations. Since natural emission Is Incoherent the third and fourth 
Stokes parameters are zero. It Is usual to assume that T* are Independent of 
the azimuth angle CIO]. Thus, 
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where P(±u s »u) Is the zeroth-order Fourier component of the phase matrix; K a 
Is the absorption coefficient matrix (Appendix A). The boundary condition of 
the ground surface Is 


T + H,u s ) ( G(u s , 
Jn 


u) T“(-d,u) dy + e g T g , 
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where 6 Is the zeroth-order Fourier coefficient of the surface phase matrix, 
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; s is the emlsslvlty of the ground surface, end T g ,s the temperature of the 
ground medium The procedure to solve (.)-(!*, up to the first-order , spuming 

7 5 S “ e " ,n the !«*'- 9lv.n in the previous section 
when . leave out step The solution for T- ( o,) including seroth-order 
and significant terms In the first-order Is 
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Where 1 Is the Identity matrix; . is the albedo of the layer; are the 

outgoing and Incoming direction cosines corresponding to the points chosen for 
an h-polnt quadrature Integration; and e’ e1 d 1s the diag0 „ a , 

- Dlag [^1 e -K eh d/ Ul . 

In (13), Hi, j<N. The matrices M 1(2 3 represent volume 
and their elements denoted by (^) for a polarisation component p, (p.,., 
or h) are related to the elements o/the f , s 
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scatterl ng operators 


C" 2 p q ) 


'ij * 0,5 K sp1 f p pq tu 1 '“j 5 + P pq(“i>-'‘j)) (l - e" Ke P ,d )/K 


epi 


r 

i 

i 


(%))„ * °- 5 K spi' p p,<“l-4) t 1 - e ’ (KePl+Keqj)d ] / ( K epf%j) • 


I. where K sp1 »K sp /u 1t K e p1“ K ep /u 1» and K eq j»K eq /uj. In (13) terms Involving M 
j represent the first-order corrections to the zeroth-order terms. It Is 

1 Interesting to note that Irrespective of order the major sources of emission 

| are the upward and downward emissions of the layer and upward emission from 

the ground. These source emissions can either propagate directly or be scat- 
I tered by the layer boundary towards the point of observation with an appro- 

| priate amount of attenuation (zeroth-order terms). They can also be scattered 

by the layer Inhomogeneities with or without additional scattering by boundary 
j towards the point of observation (first-order terms). Terms which represent 
scattering of these emissions first downward and then upward by the ground 
i surface have been ignored In (13). 

| Usually the upward emission from the layer Is at least comparable to that 

from the ground for a fully grown vegetation cover. Hence, the ground 
| emission term cannot dominate the total emission. The downward layer emission 

which Is scattered upward by the ground Is also sensitive to soil moisture 

» condition. Its Increase with the Increase In soil moisture compensates the 

% * 

J simultaneous decrease In ground emission. As a result the sensitivity to soil 

moisture can be reduced by 50% or more when the optical depth of the layer Is 
| 0.4. 
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4.0 SCATTERING CHARACTERISTICS OF A VEGETATION LAYER 

In this section we want to show the relative contributions of the terms 
in (5) as a function of the incidence angle at different albedo, optical 
depths, and ground permittivities. For the purpose of Illustration, Klrchhoff 
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surface model written In series form [11] and Rayleigh phase matrix [2] are 
used In Figs. 2-5. Fig. 2, Illustrates the contributions of the three types of 
terms discussed In Section 2 at an optical depth of 0.1. It Is seen that the 
surface-volume scattering term Is much more Important for horizontal 
polarization than vertical polarization at large angles of Incidence. Fur- 
thermore, the horizontally polarized surface-volume scattering term also 
becomes comparable to the volume scattering term at large angles of Inci- 
dence. This Is true particularly when the ground moisture Is large while 
optical thickness Is relatively small. In this case, soil moisture sensing Is 
possible at large angles of Incidence even though the surface scattering term 
Is small at these angles. Fig. 3 shows that for thicker optical depth the 
surf ace- volume scattering term becomes comparable to the surface scattering 
term at much smaller Incidence angles than when the optical depth Is small 
(Fig. 2). This Implies that the surface-volume scattering term tends to 
reduce the loss In sensitivity to soil moisture sensing as optical depth 
Increases. In Fig. 4 the backscatterlng coefficients for the three types of 
terms In (5) are plotted versus the optical depth at two different permit- 
tivities when the Incidence angle is 8.6°. The surface-volume term Is small 
compared with the surface term. When the same calculations are repeated at 
20° Incidence In Fig. 5, we see that the surf ace- volume scattering term 
becomes Important and can exceed the surface scattering term If optical depth 
Is large. Since soil moisture sensing Is usually conducted between 10° to 
20°, these results Indicate that surface-volume scattering terms tend to 
reduce the loss In sensitivity of nhe surface scattering term when the optical 
•depth Is large (t>0.4). To summarize these results, the ratio of the change 
In the backscatterlng coefficients between a wet (soil permittivity * 18) and 
a dry (soil permittivity * 3.5) soil condition, with vegetation cover present 
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to that with no vegetation cover Is plotted In Fig. 6 using the disc phase 
function (Appendix A). Th,1s ratio Is shown as a function of optical depth. 

It depicts the loss In sensitivity In sensing a vegetated soil versus a bare 
soil. As expected the loss In sensitivity Increases with the Increase In 
optical depth or albedo. 

5.0 EMISSION CHARACTERISTICS OF A VEGETATION LAYER 

The angular characteristics of emission from a wet (c«25) and a dry soil 
(e«4) surface with vegetation cover is shown In Fig. 7. Unlike surface emis- 
sion, there Is very little difference In angular trends between horizontal and 
vertical polarizations at an optical depth close to unity. Note that the 
zeroth-order theory [(13) without the -terms] In common use Is quite Inade- 
quate unless both albedo and optical thickness are very small. This point Is 
Illustrated In Fig. 8 where zeroth, first, and exact (numerical) solutions are 
compared. 

In Fig. 7 the level of upward emission by the layer Is much higher than 
the emissions by the soil surface at all nadir angles even for the dry soil 
case. At smaller optical depth (t<0.8) (Fig. 9) surface emission may be com- 
parable or higher than layer emission. The surface emission term Is higher 
and the I term (downward emission -scattered upward by the soil surface) Is ' 
lower for the dry case than the wet case. This shows that the change In soil 
emission due to change In moisture Is reduced by the term In general. As 
albedo decreases It Is expected that layer emission will Increase while sur- 
face emission may suffer a slight decrease. 

Thus, for the same optical depth, a smaller albedo leads to a smaller 
difference In total emission between the wet and the dry soil conditions 
because the upward layer emission Is larger and surface emission Is somewhat 
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lower. In summary, the ratio of the change In emissions from a wet and a dry 
soil surface with vegetation cover to that from the same surface without 
vegetation cover Is shown as a function of optical depth In Fig. 10 using the 
disc phase function. It Is seen that a larger albedo results In a better 
sensitivity to soil moisture change because scattered contribution from soil 
surface has Increased while the upward layer emission has decreased. Upon 
comparing Fig. 6 with Fig* 10, It Is seen that, In general, vegetation causes 
a more severe degradation In sensitivity on passive sensing than on active 
sensing. It Is Interesting to note that while larger layer albedo helps to 
Improve the sensitivity somewhat In passive sensing, It has a significant 
adverse effect on active sensing. 

Summary of Layer Parameter Effects On 
Soil Moisture Sensing 


Layer 

Sensitivity 

Sensitivity 

Parameters 

Active 

Passive 

. / 

N 

/ 

T / 




6.0 COMPARISON WITH MEASUREMENTS 

The reductions In sensitivity for the active (Fig. 6) and passive (Fig. 
10) sensing of soil moisture are computed from reported measurements and 
plotted In Figs. 6 and 10, respectively. The procedure used to generate the 
data points In Fig. 6 Is as follows: 

(a) Scattering coefficient, o°, measurements at 10° Incidence angle, 
C-band, HH polarization, on bare soil, wheat, corn and soybeans 
are taken from [12,13,14]. o° measured from wheat, corn, and 
soybeans will be referred to as o° (total). 
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(b) Regression lines are generated for each crop and bare soil In the 
form of o° versus soil moisture. An example Is shown In Fig. 11. 

(c) To estimate optical depth at 10° Incidence for a given crop, we 
plot o°(total) versus o°(bare soil) for various values of soil 
moisture [6]. This plot Is In real numbers. Since o°(total) a 
[exp(-2x/cos0)] o°(bare soil) + o°(crop), the slope of the re- 
gression line of o°(tota!) versus o°(bare soil) Is equal to 
exp(-2x/cosfl). Hence, the optical depth x can be estimated [6]. 

(d) Knowing the optical depth, we can use the regression lines for 
o°(total) and o°(bare soil) to find a value for the ratio, 

["Vwrtl -J°dBH tota1 

t’V"*) - ° 0 d B H bare I011 

which Is plotted as a data point on Fig. 6. 

The data points In Fig. 10 are obtained from the emisslvlty data reported 
by Wang et al. [15]. The procedure Is as follows: 

(a) The regression lines on the emlsslvity data obtained by Wang et 
al. [15] for various vegetations versus moisture (or permit- 
tivity) are fitted by the theoretical model [(13)]. It turns out 
that there is no appreciable difference between the theoretical 
and the regression lines. Hence, each line In Fig. 12 represents 
both regression and theoretical lines. It Is found that the 

• slope of the regression line Is sensitive to the optical depth, 
while the level 1$ sensitive to albedo. Hence, each fit provides 
an estimate of both albedo and optical depth (Fig. 12) for a 
given vegetation layer. 

(b) Knowing albedo and optical depth for a given vegetation layer we 


13 




can compute the ratio 

[T(w«Q - T(dry)] total 
[T(v»t) - T(dry) ) bare J0t) 

and plot It as a data point In Fig. 10. 


7.0 CONCLUSIONS 

Scattering and emission theories up to the first-order are derived and 
expressed In a form which permits easy Identification of the zeroth- and 
first-order terms as well as simple Interpretation of the scattering and 
emission processes. Upon comparing the loss In sensitivity In soil moisture 
sensing due to vegetation cover, It Is found that the loss Is less In active 
sensing than passive sensing. One reason Is that the volume-surface 
scattering term contributes positively In active sensing, while the downward 
emission scattered upward by the ground surface contributes negatively In 
passive sensing. Another reason Is that emission from vegetation Is of the 
same order as that from the ground while scattering from the ground Is usually 
much larger than that of the vegetation at small Incidence angles. 
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APPENDIX A: Disc Phase Function and Extinction Coefficients 




The phase matrix for a dielectric thin circular disc may be written In 
terms of the scattering amplitude matrix f [16] given as 


where 


f ■ [k 2 a t (c-1) Oj (a o t )/(2 q t )] 


* A T A 

(A*v)‘*v s 

m _ 

A T A 

(A* v) »h s 


* * T * Tl 

(A»h) *v 

d 

m _ 

A T A 

(A*h) 1 *h s 
(A-J ) 


k * 2ir/\ : free space wave number, 
a : radius of disc, 
t : thickness of disc, 
e : permittivity of disc, 

: first-order Bessel function of the first kind, 
q t «k[(s1ne s cos$ s -s1n9cos<f>r+(s1ne s s1n$ s -s1nes1n<fr) ], 
e $ ai:d 0 are polar scattered and incident angles, respect 1 vely, 

♦ s and $ are azimuthal scattered and Incident angles, respectively, 
(see Fig. 1 for geometry) 

AAA A 

v, h, v s , and h $ are unit polarization vectors as shown in Fig. 1, 

A 

h « (-sin*, cos<fr, 0) 
v ■ (cosecos*, cos0s1n$, -sine) 
h s » (-s1n* s , cos* s , 0) 

A 

v $ « (cose $ cos<fr s , cose s s1n<j> s , -s1ne $ ) 

* 

A Is a three by three matrix given as 




A » ■— 


.0 0 a ()/ a lJ 


where 


, 2 . 

a n ■ 1 + V- («-l) A f 


■ 1 + ^ 2 ^ (e-1) Aj 


(A-2) 


A 


0 


A 


1 


(a 2 -c 2 ) 


-1.5 


2(a 2 -c 2 ) 


-1.5 




0.5 


LaVl 


0.5 



In the above, a and c are major and minor semi -axes of an ellipsoid. Note 
that the symbol T In (A-l) denotes transpose of a column matrix. 

From (A-l), the phase matrix P may be constructed as 

P * 4 ir nQ <o> 

where o Is the Stokes matrix expressible In terms of f [9, p« 35] and n Q Is 
the number density of discs. <> denotes averaging over three modified 
t'ulerlan angles [17] which account for the orientations of a disc. 

The extinction coefficient matrix K e Is given as 

* diag(K v , K^, K^, K^) 
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where 


K h * k ah + K sh 
K 3 • K 4 - 0.5(K V +K h ) 

The scattering coefficients K sv and are 

K sv ■ n 0 J di! < l f vvl 2 + l f hvl 2> 

K sh * "0 ] d!l <l f hhl 2 + l f vhl 2> 

where / dn denotes Integration over the solid angle 4 tt. The absorption coef- 
ficients K av and are [9, p. 17] 

K av * n o k E " * ^ 1 < l A#v -jl^ > 

K ah * n 0 k e " * 1 <|A*h^| 2 > 


where e" Is the Imaginary part of the disc permittivity c 


FIGURE LEGEND 


Figure 1. Geometry of the scatter and emission problem. 

Figure 2. Angular behavior of backscatterlng coefficient of Rayleigh layer 
above Irregular ground. (w*0.1, x*0.1, ko*0.1, k£*10; where <o : 
layer albedo, x : layer optical depth, o : standard deviation of 
rough surface height, £ : correlation length of rough surface 
height, e g : ground permittivity, and T * S+I+V, where S : surface 
scattering term, I : surface-volume scattering term, V : volume 
scattering term). 

Figure 3. Angular behavior of backscatterlng coefficient cf Rayleigh layer 
above Irregular ground. (w*0.36, x*0.98, ko*0.1, k£*10; and T * 
S+I+V, where S : surface scattering term, I : surface-volume 
scattering term, V : volume scattering term). 

Figure 4. The effect of optical depth on HH-polarized backscattering 

coefficient at the Incidence angle 0*8.6°. (w=0.3, ko=0.5, k£*10; 

and T * S+I+V, where S : surface scattering term, I : surface- 
volume scattering term, V : volume scattering term). 

Figure 5. The effect of optical depth on HH-polarized backscattering 

coefficient at the Incidence angle 0*20°. («*0.3, ko*0.5, k£*10; 

and T * S+I+V, where S : surface scattering term, I : surface- 
volume scattering term, V : volume scattering term). 


Figure 6. 


Figure 7. 


Figure 8. 


Figure 9. 


Sensitivity reduction In backscatterlng due to a vegetation cover 
of optical depth t. (ko*0.5, kz*10; and Sensitivity Reduction * 

[®dB <wet) ' 0 dB (dry) ltotal / I°dB (wet) ' °dB (dry) lbare soil' wet 
ground c«18; dry ground e*3.5.) The disc phase function Is used: 

disc thickness * 0.2 mm, disc permittivity * 23.4+j4.7, frequency 

* 4.3 GHz; case (1) radius of disc * 1.3 cm, volume fraction of 

discs • 0.35% (equivalent albedo ■ 0.3); case (11) radius of disc 

■ 0.75 cm, volume fraction of discs « 0.2% (equivalent albedo * 

0 . 1 ). 

Angular behavior of emission of Rayleigh layer above Irregular 
ground. (u*0.36, t« 0.98, ground rms slope * 0.15, ground (soil) 
and Rayleigh layer temperature * 290°K, and T*V+I+S, where V : 
layer upward emission, I: layer downward emission scattered upward 
by soil, S : soil emission scattered and attenuated by layer). 
Comparisons between numerical, first-order, and zeroth-order 
solutions In emission for different albedo w and optical depth 
t. (The ground Is assumed to be a plane. Ground and Rayleigh 
layer temperature * 300°K.) 

The effect of optical depth on brightness temperature at the nadir 
angle 0*5.6°. (w*0.3,' and T*V+I+S, where V : layer upward 

emission, I : layer downward emission scattered upward by soil, S 
: soil emission scattered and attenuated by layer, ground and 
Rayleigh layer temperature « 290°K.) 
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Figure 10. Sensitivity reduction In emission due to a vegetation cover of 

optical depth t. (« : albedo, 0 : nadir angle, ground rms slope * 
0.15, Sensitivity Reduction « 

[T(wet) - T(dry)] tota ,/[T(*et) - T(dry)] bare so1 ,, wet ground 
e»18, dry ground e«3.5, T : horizontal polarized brightness 
temperature. ) The disc phase function with the same parameters 
as In Fig. 6 Is used. 

Figure 11. Backscatterlng coefficient as a function of soil moisture content 
mf for wheat field at C-band, 10° Incidence angle. 

Figure 12. Linear regression fit of emissivlty from different vegetation as a 
function of soil moisture content m y . (Albedo w and optical depth 
t are estimated by fitting the theory to a given regression fit. 
The ground permittivity e»3.5 and 18 are assumed to correspond to 
the soil moisture m v «7% and 21%, respectively.) 
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